from scipy.integrate import solve_ivp

# Define the ODE system
def system(t, y):
    return [y[1], -y[0]]

# Solve the ODE
sol = solve_ivp(system, [0, 10], [1, 0])
print("Solution:", sol.y)